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Abstract. We investigate the ground state of a one-dimensional lattice system that 
hosts two different kinds of excitations (species) which interact with a power-law 
potential. Interactions are only present between excitations of the same kind and 
the interaction strength can be species-dependent. For the case in which only one 
excitation is permitted per site we derive a prescription for determining the ground 
state configuration as a function of the filling fractions of the two species. We show 
that depending on the filling fractions compatible or incompatible phases emerge. 
Furthermore, we discuss in detail the case in which one species is strongly and the 
other one weakly interacting. In this case the configuration of the strongly interacting 
(strong) species can be considered frozen and forms an effective inhomogeneous lattice 
for the other (weak) species. In this limit we work out in detail the microscopic ground 
state configuration and show that by varying the density of the weak species a series of 
compatible-incompatible transitions occurs. Finally we determine the stability regions 
of the weak species in the compatible phase and compare it with numerical simulations. 
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Understanding the ground state properties of one-dimensional lattice gases or spin- 
models with long-range interaction is a notoriously difficult problem in many-body 
physics. Only few analytical solutions for classical nua and quantum j31 H] El |B] models 
have so far been found. Moreover, quantum systems with long-range interactions are 
typically difficult to treat numerically even with modern numerical techniques based on 
the density-matrix renormalization group EH Ei m no]. 

Despite these difficulties physical systems with power-law interactions have received 
significant attention in the last decades. In early years the the Kondo problem for 
the one-dimensional Ising model with inverse square interactions m, and the phase 
diagram thereof [T5] were analyzed with renormalization group techniques. In more 
recent years the critical properties of spin models with general algebraic interactions 
and dimensionality were studied mmm- 

This work connects to the studies in Refs. DEI where a classical one-dimensional 
Ising model with an external magnetic field and convex (long-ranged) interactions was 
investigated, for which the ground state was constructed analytically. In Ref. [IB] it was 
moreover shown that here the ’spin-up’ excitations, upon varying their density, form a 
complete devil’s staircase. The role of quantum fluctuations in these spin models was 
explored for the case of power-law potentials of the form l/r a in the context of dipolar 
Bose gases HZ! (a = 3) and strongly interacting Rydberg gases [18] (a = 6). Recently, 
the onset of a devil’s staircase has indeed been identified in a one-dimensional lattice 
gas of atoms that were laser-excited to electronically high-lying Rydberg states [19]. 

In this paper we generalize the aforementioned classical studies to a two-component 
lattice system. Specifically, we consider a situation in which there are two different kinds 
of excitations to which we will refer as s (strong) and w (weak). Power-law interactions 
among s(u>)-excitations are assumed to be strong (weak) while there is negligible 
interaction between excitations belonging to different components. One motivation for 
conducting this study is to show that the physics of this system is much richer than 
the single-component case and in fact features a series of transitions that are absent 
in the single-component situation. A second motivation is derived from the fact that 
multi-component lattice gases with power-law interactions move more and more into the 
focus of theoretical and experimental studies on cold atomic gases in which atoms are 
excited to Rydberg states [20; ill. j22. 23,123J [25] 26, 27]. In certain regimes aspects of the 
physics of this atomic system — specifically, strongly state-dependent interactions and 
highly suppressed interactions between atoms in different states [27] [2B] - are indeed 
well described by the model discussed here. 

The paper is organized as follows. In Section [2] we define the model. We then 
summarize the results obtained in [Q [2, [TB] in Section [3] for the single-species case. We 
do so because we will make an extensive use of the technology and formalism developed 
in those works. In Section [4] we introduce the concept of compatible and incompatible 
ground states and in Section [5] we introduce an algorithm for finding the microscopic 
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arrangement of excitations in the ground state. Finally in Section [6] we identify the 
locations of the compatible-incompatible transitions in the special case of filling fractions 
of the form 1/q, q € N + . 


2. The model 


We focus on a one-dimensional lattice system where each site is occupied by a three- 
level system whose internal states belong to the set {|0), |s), |w)}. The three levels 
correspond to the (empty) state 0 and the cases in which the site is occupied by an 
s- or w-excitation. Excitations of the same kind interact with an inverse power-law 
potential with exponent a and strength (with /a = s,w and > 0), while there is 
no interaction between species of different type. The Hamiltonian then reads 



-K n k 


l>k 


n k n l 

o i - ky 


(i) 


where n ^ = /i) (/i and h M is a positive parameter. Linking this model to experiments 
with gases of Rydberg atoms, one would think of |0) as being the atomic ground state. 
The state /i) then corresponds to an electronically high-lying Rydberg i/S-state with 
principal quantum number |29j which is excited with a laser of detuning /z /( . In 
order to have one strongly and one weakly interacting species one needs to require 
u 8 In the following, unless specified otherwise, we consider the asymptotic 

scenario V s /V w —> oo. 

Our goal is to find the configuration of the two species that minimizes the 
energy of Hamiltonian ([Tj) given fixed values of the individual total excitation numbers 
= =1 njj. . The condition V s /V w —* oo renders the search for the ground state 

into a two-step process: First, the ground state configuration of the s-excitations is to 
be found which subsequently remains “frozen 1 '. This leads to an effective — in general 
inhomogeneous — lattice on which excitations of the w-species are to be arranged in the 
second step. Depending on the imposed total excitation numbers one can distinguish 
dispositions of the w-excitations that are compatible with the one of the s-excitations, 
and others that are incompatible. Transitions between these two cases can then be 
studied as a function of the total excitation numbers n M (see e.g. [31] for a related 
discussion of compatible-incompatible transitions). 


3. The ground state of the single-species problem 

The single-species case can be described in terms of the Hamiltonian (JTj) once we set 
e.g. V w — h w — 0. Our aim is to End the configuration that minimizes the interaction 

$ In fact the strength of among excitations of the same species (atoms in the same Rydberg 
state) scales with the eleventh power of the principal quantum number Moreover, as shown by 
a perturbative calculation J2B] and recently also experimentally m there are states for which the 
interactions among excitation of different species is strongly suppressed. 
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Figure 1 . The single species devil’s staircase for a potential V(r) = 1/r 2 . 


energy once a density of excitations p = m/N is giver{§] For a start let us consider 
three excitations arranged on a line with distances r (between the first and the second 
one) and r' (between the second and the third one) with r < r'. The solution of the 
energy minimization problem is obtained by realizing that for a general convex potential 
V (r) the variation of the potential energy under a change of the position of the second 
excitation by 5r obeys 

V (r + Sr) + V{r' - Sr) < V(r) + V{r'). (2) 


This expression holds for r, Sr < r'. Hence moving the excitation in the middle towards 
the third excitation leads to a lowering of the potential energy until r = r'. The 
energetically most favourable configuration is therefore assumed when both distances 
are equal, i.e. the excitations are arranged in the most uniform distribution. This 
consideration can be extended to larger systems with excitation density p where the 
ground state is a regular arrangement of excitations with inter-excitation distance 1/p. 

On a chain, however, where the coordinates can only be multiples of the lattice 
spacing the situation is much less trivial. Here the construction of the ground 
state configuration proceeds as follows (Note, that the lattice spacing is set to one, 
throughout.): We define Xk as the position of the fc-th excitation and denote the distance 
between the k -th excitation and its nearest excitation by r{ — Xk+i — Xk, the distance 
between the fc-th and the next-to-nearest excitation by r^ = Xk +2 — Xk = r[ k) + r[ k+1 \ 
and generally the distance to the /-th-nearest excitation by r\ k ' > = = 

(k) J 

Xk+i — Xk- All spacings of neighboring excitations, r) , will belong to the set <Si = 
{[1/pJ, [1/p]} which is referred to as the minimal set. 

Using the above abbreviations the task is now to minimize 


E 


m 

—mh + 

i=i 


5>(U>) 


( 3 ) 


§ In this section we drop the label s as we consider a single species. 
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In Ref. [Ij it was proposed to look for a solution which minimizes each inner sum 
separately, knowing that if such solution exists it will also minimize E. There it was 
demonstrated that this solution exists and that it is achieved when r\ ; form minimal 
sets, for each l. These sets can be defined as Si = {ri,ri + 1}, where ri = |_ l/p \, and 

ri + 1 = \l/p\- 

Another challenge is finding the microscopic arrangements of the excitations on an 
infinite chain with rational filling fraction p — p/q. Clearly for all the filling fractions 
of the form p — 1, q G N + one just has to distribute the excitations uniformly, one 
every q sites. For general rational numbers it is more complicated to figure out how the 
excitations are distributed, but we can still think that the configuration will be periodic 
in q, each period having p excitations arranged in the most uniform way. The algorithm 
for finding such a distribution is given in Refs. mca. To illustrate it we start with an 
example and then report the general solution. 


3.1. Example 


Let us consider the case p = 11/47. From the previous explanation we know that 
the distances between neighboring excitations are r[ k ^ G {4, 5}. This automatically 
minimizes ^2 k V(r ^). If we just had to minimize this contribution we could distribute 
the 11 excitations over the 47 sites in any pattern, which would ensure that excitations 
are four or five sites apart. However, as we have to minimize each contribution l in (J3]) 
separately, we need a more sophisticated way to understand the configuration. We start 
by writing 


1 , 3 

_ — 4 + — , 
P 11 


( 4 ) 


hence finding r\ k) G Si = {4, 5}. With 47 sites we therefore find 8 intervals r of 4 sites, 
and three of 5 sites. The only possible values for the distances between next-neighboring 
excitations r 2 , are then given by 8, 9, 10. Note, that for some part of the following 
discussion it turns out that it is actually convenient to express these distances as strings 
of consecutive spacings, i.e. 44, 45, 55. Let us now find the permitted next-nearest- 
neighbor distances which form the set S 2 - This set must contain only two consecutive 
numbers, so either 8, 9 or 9,10. In order to understand which ones to pick and how to 
distribute intervals of length 4 and 5 we define a new filling fraction using the non-integer 
remainder of Q: 


— = 4 - -. 
3 3 


( 5 ) 


This filling fraction represents the density of 5 sites intervals which can be regarded 
as defects within the 4 sites intervals. We thus have two strings of length 4, and one 
of length 3 between these defects (4+4+3=11), where the strings of length 4 and 3 
correspond to 4445 and 445, respectively. 

In order to understand how these strings are distributed we iterate the logic and 
obtain from the remainder of ^ the filling fraction (1/3) -1 = 3. This means that if we 
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identify the string 445 as the defect, and 4445 as the defect-free one we have eventually 
44454445445. This is the ground state configuration as the iterative procedure stops 
when the effective 1/p G N + (here (1/3) -1 G N + ). This method creates the configuration 
which minimizes each contribution in (|3| separately. In fact one can check that for 
V/ < m , ri G Si. 


3.2. General algorithm 


The generalization of this method was given in Ref. |lj. It starts by considering the set 
of numbers 


1 

P 

1 

r 0 

1 

T 1 


= n + r 0 


= 7ii + r i 


= n 2 + r 2 


1 


T'k-l 

where Uj G Z and A: is a finite number since p G Q, and —1/2 < Tj <1/2 for any 
j. Then one can define iteratively the strings of non-defective sites X], X 2 ,..., X k and 
defects Yi, Y 2 ,..., Y k in the following way 

X-! = n, Yi — n + sgn (r 0 ) 


Tk-2 

1 


k -1 ' 1 fe -1 




X i+1 = [.X i]" 4 " 1 Y i: Y l+1 = Y t . (7) 

Here the notation [. X] n means that the spacing X is repeated n times. The ground state 
configuration is then given by X k . 

3.3. Stability of the ground state 

Finally, we summarize the solution to the problem of understanding in which region 
of the parameter space — and in particular for which values of the parameter h — a 
given configuration of density p = m/N is stable. The solution was given by Bak in 
Ref. [16j assuming a chain of length N with periodic boundary conditions. The essence 
is to analyze the energy cost of inserting or extracting one excitation from a given 
configuration. If both operations increase the energy (J3| the configuration is stable. 

Adding one excitation will modify the number of Z-tli nearest neighbor distances r; 
and T\ T 1. Starting from a given configuration with m — a distances r/ and a distances 
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T'i + 1 the following system of equations must hold 

f EZi r i k) = ( m - a )c + «(n + 1 ) = ni 

\ E7Y\ r< t k) = (m + 1 - a')ri + a'(n + 1) = iVZ 

Solving ([ 8 ]) gives the number a' of distances r t + 1 in the configuration with one added 
excitation, a' = a — ri. It means that 77 distances ri + 1 will be replaced by ri + 1 
distances 77 . Note that in the special case mri = Nl , one instead has to replace 77 — 1 
distances 77 by ri distances 77 — 1. With the same logic one understands that extracting 
one excitation from the system will replace 77 + 1 distances 77 with 77 distances 77 + 1. 
Considering the energy ([ 3 ]) we have 

A E + = —h+ + (ri + l)y(n) - t’iV(t’i + 1) + (r 2 + l)V(r 2 ) - r 2 V(r 2 + 1) + ... + 

+ qV(q — l) - (q— 1 )V(q) + ... + 2 qV (2 q - 1) - (2 q - l)V(2q) + ... 

A£L = h_ - (n + l)C(rx) + riV(r x + 1) - (r 2 + l)C(r 2 ) + r 2 V(r 2 + 1) + ... + 

- (q + 1 )V(q) + qV(q + 1) + ... - (2g + 1)V{2q) + 2gC(2g + 1) + ... 

(9) 

where we are denoting with A E± the change in energy when adding or extracting an 
excitation. Setting the two quantities A E± equal to zero gives a region of stability of 
the width 

h + — h_ = ^ nq [V(nq — 1) + V ( 7 iq + 1) — 2V(nq)\ . (10) 

n 

In [TBj it was proven that these intervals of stability for rational values of p fill up the 
whole parameter space h G R + , giving rise to a so called devil’s staircase as shown in 

Fig{l| 

4. The two-species problem — Compatible and incompatible ground state 
dispositions 

We will now turn to the case in which there are two different species — s and w - 
in the lattice. If we allowed for multiple occupations on each site we would clearly 
have that the energy (jl| is minimized by disposing s and w excitations independently, 
using the algorithms (|6]) and ([7]). However, as we allow just for one excitation per site 
we have the additional restriction p s + p w < 1. Nevertheless, we will see, that in some 
cases, i.e. for some filling fractions, the two populations do not influence each other. The 
more interesting case, however, is encountered when both species cannot simultaneously 
assume their ground state configuration. We will refer to the former case as compatible 
and to the latter as incompatible. 

The goal of this section is to define a method for distinguishing the compatible 
cases from the incompatible ones. Note that this is a purely combinatorial problem and 
that energy considerations are irrelevant for the notion of compatibility. To distinguish 
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Figure 2. Compatible (black squares) and incompatible (white squares) filling 
fractions for a periodic lattice of length N = 16, 20,47. The results are given for 
different excitation numbers n s and n w corresponding to filling fractions p s = n s /N 
and p w = n w /N. The shaded area corresponds to cases p w + p s > 1 which cannot 
occur since each lattice site only contains at most a single excitation. 


compatible from incompatible cases we consider the filling fractions p s = n s /N and 
Pw = n w /N, such that the distributions of the s- and w-species are periodic with the 
same period N. This causes no loss of generality as the filling fractions with two different 
periods, e.g. p s = n s /N s and p w = n w /N w , can be always expressed in terms of a 
single period, p s = n s N w /(N s N w ) and p w = n w N s /(N s N w ). The two configurations 
corresponding to p s and p w are compatible if within a period each excitation of the kind 
w is located in a hole with respect to the arrangement of the s-excitations (s-holes). The 
density of s-holes, denoted by p s , is linked to the density of s-excitations by p s = 1 — p s , 
and the disposition of such holes can be evaluated by using the algorithm (Jt|) . As in 
the previous section we start with an example and then provide the general rule for 
distinguishing a compatible case from an incompatible one. 

4-1. Example 

Consider the filling fractions p w = 3/16 and p s = 9/16. The latter yields p s = 7/16. 
The disposition of w-excitations is 556, while the disposition of the s-holes is 3223222. 
If one is able to match the positions of the 3 ic-excitations with the positions of three 
of the 7 s-holes without distorting the excitation distance pattern 556, then the w and 
s configurations are compatible. One can therefore think of the hole configuration as 
an effective inhomogeneous or distorted lattice with spacings a = {2,3}, disposed with 
repeated periods ...PPP..., P = 3223222. Using ([7]) with the effective filling fraction 
Pw/ps = 3/7, one finds that the configuration of the 3 w- excitations on the distorted 
lattice is 223. Here the first 2 in the string 223 means the first two elements in P 
and so on, i.e. disposing the w-excitations on the distorted lattice gives the disposition 
223 —> (32) (23) (222) = 556. This is identical to the original disposition, so that the 
dispositions of the w- and s-species are compatible. With the same logic one finds 
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that the filling fractions p w = 3/16, and p s = 12/16 are not compatible. In this case 
p s = 4/16, with a configuration of holes 4444, and there is no way we can combine 
multiples of 4 in a way to match the distances 5 and 6 between ic-excitations. 


4-2. General algorithm 


The general algorithm to check the compatibility of two distributions of excitations is 
based on the following sets of numbers 


Vi = {r h ri + 1} 

Ui = {s h si + 1 } 

Qi = Ui,qi + 1 } 



(ii) 


These are the distances Vi between /-th-nearest w-excitations and the distances T~Li 
between the /-th-nearest s-holes. The third set contains the distances between /-th- 
nearest u’-excitations counted in units of lattice spacings of the distorted lattice. We 
focus here on the case p w /p s < 1 (i.e. p s /p w > !)• The opposite case is described by 
simply inverting the roles of excitations and holes. 

The distance between Z-th nearest w excitations is either qi or qi + i counted in 
distorted lattice spacings, where a distorted lattice spacing reads either si or si + 1 
original lattice spacings. If it is possible to group qi and qi + 1 distorted lattice spacings 
Si and Si + 1 in a way that they sum up to give the distances rj and r; + 1 then the 
two densities are compatible. For example, the compatible case in Fig. [3] has p s — 1/2 
and p w = 1/4. Notice, that in this simple example the sets Eq. GD contain only one 
number and read Vi = {4/}, Qi = {21} and for / = 1 "Hi = {2}. Clearly, the qi = 21 
distances s\ — 2 yield the distance 4/ = ri, i.e. the two densities are compatible. On the 
other hand, considering e.g. the first incompatible filling fractions in Fig. [3] = 5/16, 

p w = 1/4 one has for / = 1 V\ = {4}, "Hi = {1,2}, Q\ = {2,3}, and in particular the 
distorted lattice is formed by repeated periods P = 12121212112. One can group 3 
consecutive distances in P to get r\ —4 (namely 112), but it is not possible to obtain 4 
from grouping only 2 consecutive distances, so that already at level / = 1 one can prove 
the incompatibility. Formally, this means that two configurations are compatible if 


VI < n w , H qi O / 0 A P qi +i nVi^Q. (12) 

In Fig|2] we show a diagram illustrating the filling fractions that lead to compatible or 
incompatible dispositions of excitations for different system sizes. 

As a byproduct of Eq. (12) one learns that a sufficient condition for two filling 
fractions to be compatible is that they can be represented as p w = n w /N, and 
p s = 1 — kp w , with 1 < k < [p~ l \. In fact when k — 1 the two configurations are dual 
(corresponding to the upper left to lower right diagonal in Fig. [2]), in the sense that 
Pw — Ps- Increasing k increases the number of holes while the w- and s-configurations 
remain compatible. 
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Figure 3. (Color online) Microscopic configuration for different choices of the Riling 
fractions p s and p w . The excitations of the s-species are indicated with solid blue 
large circles and those of the ir-species with solid red small circles. Empty sites are 
represented by empty circles. The case shown at the bottom of the figure is compatible 
(green shade), such that the arrangement of both the s- and w-excitations can be 
achieved by applying the algorithm 0 on the two species independently. All other 
examples are incompatible (red shade). 


5. Disposition of the excitations in the incompatible case 

Let us now discuss the actual ground state of the Hamiltonian (JTj) . For compatible filling 
fractions this is given by the distributions of the two species arranged independently 
using the algorithm (JT]) . In case of the two filling fractions being incompatible the 
situation is more complicated. In order to make analytical progress we work under 
the assumption that V s /V w —> oo, suchthat the arrangement of s-excitations can be 
considered as “frozen”. Note, that such situation can be achieved in the context of a 
realisation of the system with atomic Rydberg gases when the two principal quantum 
numbers corresponding to the s, w states are chosen appropriately, see for more 
details. The problem then reduces to the arrangement of ic-excitations interacting with 
a l/r“ potential constrained to sit in the s-holes. It is now convenient to define an 
effective inhomogeneous or distorted lattice with spacings a G , f 1 ^. Clearly 

the inhomogeneity alters the minimum energy disposition of the w-excitations, such 
that we cannot consider only nearest neighbor distances rq and rq + 1 , but we have in 
general to consider distances up to rq + /?, with (3 > l. In the following we use the term 
“distortion of the minimal set” to denote that Z-th nearest neighbours distances can take 
values in the set { 77,77 + 1, ...,77 + /?} rather than { 77,77 + 1}. The smallest distortion 
then refers to the case when only one more distance ri + 2, i.e. f3 = 2, is included. 

One can now exploit the convexity of the potential to demonstrate that the energy 
is minimized by the most homogeneous configuration. We start by considering the 
smallest distortion of the minimal set. Then the following relation holds: 

Tlw 

y j y = {n w - a- b)r t + a(r t + 1) + b(ri + 2) = IN, (13) 

i=1 

where a and b are the number of occurrences of the distances 77 + 1 and 77 +2 respectively. 
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We now want to show that the introduction of a longer distance results in an 
increase of the energy. To this end let us consider the set {77, ri + 1, 77 + 2 ,77 + 3} where 
only one additional distance 77 + 3 is introduced in the ground state configuration, while 
the total number of excitations n w is kept fixed. Then, since the sum (13) must be 
conserved, we have!"" 


' " UJ 

y, r \ l) = (n^ —a — 6 — 1)77 + 0(77 + 1 )+ (6 — 1) (77+ 2 ) +(77 —1) + (77+ 3 ). ( 14 ) 


i= 1 


Calling E 2 the energy of the configuration (13), and E 3 the one of (14) the following 
holds 


E 3 - E 2 = V(ri - 1) + V (77 + 3) - V ( 77 ) -V(r t + 2) > 0. (15) 

where the positivity stems from the convexity of the potential (J2|. Analogously, one 
can prove by induction that having a larger set of distances is always energetically 
unfavourable, such that the minimal energy configuration is indeed the one achieved 
with the minimal set of distances, i.e. with j3 — 2. 

The procedure for finding the ground state configuration in the incompatible case 
can now be summarized as follows. 

• Define the inhomogeneous lattice using ([7]) with p s . 

• Distribute the w-excitations in the effective lattice given by the disposition of s- 
holes, using algorithm (JT]) with p w /p s as the effective filling fraction. 

In Fig. [3] we show some examples of the ground state configurations for various filling 
fractions. 

We now want to show that the above-described procedure automatically yields the 
minimal distortion, i.e. /3 — 2, of the minimal set of distances of the w-excitations. Lets 
denote the two values assumed by the s-hole distances a as a',a' + 1. It follows from 
([7]), that the ground state configuration always contains strings of the form a ,k (a' + 1) 
and a ,k+1 (a' + 1), or (a' + l) fc a' and (a' + l) fc+1 a'. The ground state configurations of 
the s-holes can then be divided into four mutually exclusive cases, depending on what 
kind of string of distances a' they contain. If they contain 

• a' k+1 {a' + l)a /fc (a' + l)a' fc+1 (a' + 1), 

the maximum distortion occurs for the configurations with Qi — {k + 1, k + 2} for 
Ah-nearest neighbor I 11 this case s qi G {a ,k+1 ,a ,k (a + 1)} and s q;+ i G {a' k+l {a' + 
1), ( a' + 1 )a' k {a' + 1)}, such that the maximum difference in length is /3 = 2. 

• a' k (a' + l)a /fc+1 (a' + l)a'V + 1), 

the maximum distortion occurs for the configurations with Qi = {k + 1, k + 2}. 
The maximum difference in length is (3 = 2. 


|| Adding the distance ri+ 3 leads to new number of distances rj, n + 1, n + 2 which are constrained to 
sum up to Nl together with the single distance n + 3. This can be achieved by adding and subtracting 


ri + 3 to ( |13| ) and absorbing the — (r; + 3) term in the remaining terms. In (14) we have just absorbed 
the — (77 + 3) term in a way that makes the convexity of the potential (pi) manifest. 
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Figure 4. (Color online) Stability regions of the ui-excitations. Numerical results 
for a = 2 on a chain of length N = 16 are shown as blue dots. The density of s- 
excitations is fixed to p s = 1/2, and we set the parameters h s = V w = 1 and V s = 100. 
The superimposed red squares are the predictions of Eq. (19) for the compatible to 
incompatible transitions in the thermodynamic limit. 


• (a' + l) fc a'(a' + 1 ) k a'(a' + 1 ) k a' 

the maximum distortion occurs for the configurations with Qi = {k, k + 1}, giving 
again /3 = 2 as the maximum length difference. 

• (a' + l) fe+1 a'(a' + l) fc a'(a / + l) fc+1 a', 

the maximum distortion occurs for the configurations with Qi = {k,k + 1}. Also 
in this case the maximum difference in length is f3 = 2. 

It follows that apart for the compatible cases for which /3 = 1, the algorithm described 
above yields automatically the minimal distortion, that is (3 = 2, and consequently the 
resulting configuration minimizes the interaction energy. 


6. Stability region of the compatible phase and compatible to incompatible 
transitions 


In the case of a single species the convexity of the potential is sufficient to determine 
the stability of a certain density through Eq. (J8]) . This is generally not true anymore 
in the two-species case and it is very challenging to determine which configuration or 
disposition is stable for given values of h w and h s . 

We can make progress by again assuming that V s /V w —* oo. Here the results of the 
single-species problem apply for the s-excitations, while (13) holds for the w species, 
such that e.g. extracting one excitation leads to the following set of equations 

(n w - a - b)r t + a(r z + 1) + 5(r z + 2) = IN 

( n w — a! — b' — 1 )r z + a'(n + 1) + b'{ri + 2) = IN 1 ' 

The aim is to find the variables a', b 1 in terms of the variables a, b, similar to the case 
of a single species. However, the system of equations (16) does not yield a complete 
solution. Instead it provides the consistency relation 

a! — a — ri — 2{b' — b). 


(17) 
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In order to proceed we restrict our analysis to the case p s = l/q s , and p w = l/q w . 
The independent configurations are clearly homogeneous with periods q s and q w . The 
two configurations are incompatible if and only if q s and q w have no trivial common 
divisor. This can be shown as follows. Let us denote the positions of the excitations 
as x s = kq s and x w = lq w + 6, with <5 G Z. This means S is a constant shift of the 
^-configuration with respect to the s configuration. If x s ^ x w for VZ, k G Z, then the 
two configurations are compatible. We thus have 

5 ^kq s - lq w = p(kq' s - Iq'J (18) 


where in the second equality we explicitly collected p, the common divisor of q^, such 
that q^ = pq'p, p = s,w. 


The term in brackets at the r.li.s. of (18) is an integer number, such that for non 


trivial p ^ 1 one can always find a shift <5 (different from multiples of p) that yields 
non-conflicting configurations. On the other hand, if p — 1, such shift does not exist 
and the configurations are incompatible. 

Particularly interesting is to understand the transition from a compatible 
disposition (or phase) to an incompatible one. This can be studied as in Eqs. ©- 


(10) by perturbing the minimum energy configuration in the thermodynamic limit by 


subtracting or introducing one ic-excitation. If we perturb a compatible configuration 
(where q SjW have a common divisor p) this will with certainty lead to an incompatible 
disposition, as it will introduce q w distances q w ± 1: Every q w ± 1 distance shifts the 
excitations to the left/right by one site with respect to the ground state configuration 
(which contains only distances q w between nearest neighbours). This occurs q w times 


and since q w > p, the condition ( 18 ) cannot be satisfied. 

ls G N + for which the transition 


Let us focus now in detail on the case v = q w /q s 
points can be found analytically. Using the procedure described in Section [5] one finds 
that the set of Z-th nearest neighbour distances {lq w } is modified in the following way. 
The insertion of one excitation introduces Iv distances q w — 2 and l(q w — 2u) distances 
q w — 1. On the other hand subtracting one excitation will introduce lu distances q w + 2, 
and l(q w — 2 u) distances q w + 1. Notice that this is in agreement with the consistency 


condition dictated by (17). 

Calling h ^ the values of h w for which these compatible to incompatible transitions 
occur, we find 

OO 

hi = ^2 K 1 + - v))V(lq w ) - l(q w - 2 u)V(lq w + 1) - luV(lq w + 2)] , 

Z=1 

oo 

h ™ ] = [( X - 1 (9^ - ”)) V fe) + Kq* - 2 u)V(lq w - 1) + lvV(lq w - 2)] . 


i =l 


(19) 


Considering a general form of the filling fractions it is very hard to read a repeated 
pattern of distances, such that we have not been able to determine regions of stability 
analytically. 
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In order to verify Eqs. (19), we have determined the stability regions numerically 
for a finite chain of length N = 16, where we have fixed p s = 1/2. The results are 
shown as blue dotted curve in Fig. [4| In this particular case the only accessible Filing 
fractions are of the form m/N , m < N. However, since Eqs. (19) are only valid for filling 
fractions of the form 1 /q w , the only accessible values which can be used for a comparison 
are p w = 1/8,1/4,1/2. These regions are delimited by the red squares in Fig. [4j We 
can observe that the predictions in the thermodynamic limit tend to overestimate the 
values of the external field at which the transition happens compared to the finite chain 
calculations. This discrepancy is more pronounced for higher density of ic-excitations. 
This is understandable as higher Filing fractions introduce larger energy corrections in 
the Fnite chain. 


7. Conclusions 

We have investigated the statics of a two-component lattice gas with species-dependent 
1/r" interactions in a one-dimensional lattice. The motivation behind this study is a 
link between this system and the current experiments in which cold atoms are excited 
to multiple high-lying Rydberg states. We found that the ground state arrangement 
of the species in the lattice falls into one of two possible categories which depend on 
the Filing fractions of the two species. There is a compatible case in which the two 
species can be considered as independent, and incompatible one in which this is not the 
case. We have deFned the criteria for compatibility of the two species conFgurations. 
In the limiting case in which we consider one of the two species as frozen, we showed 
how to determine the ground state conFguration in the incompatible phase. Finally, for 
Filing fractions of the form 1 /q, we determined the stability regions of the compatible 
phases in the thermodynamic limit and compared the analytical result with Fnite size 
numerical simulation. 
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